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We use astrophysical and atomic clock tests of the stability of the fine-structure constant a, 
together with Type la supernova and Hubble parameter data, to constrain the simplest class of 
dynamical dark energy models where the same degree of freedom is assumed to provide both the 
dark energy and (through a dimensionless coupling, £, to the electromagnetic sector) the a variation. 
We show how current data tightly constrains a combination of £ and the dark energy equation of 
state wo- At the 95% confidence level and marginalizing over wo we find |£| < 5 x 1CU 6 , with 
the atomic clock tests dominating the constraints. The forthcoming generation of high-resolution 
ultra-stable spectrographs will enable significantly tighter constraints. 


I. INTRODUCTION 

The discovery of cosmic acceleration from luminosity 
distance measurements of Type la supernovas 00 led to 
wide-ranging theoretical and observational efforts trying 
to understand and characterize it. While a cosmological 
constant A remains the simplest explanation consistent 
with most observational data, the well-known fine-tuning 
problems associated with this solution imply that alter¬ 
native scenarios should be sought and tested. The most 
natural alternative would involve scalar fields, an exam¬ 
ple of which is the recently discovered Higgs field 01, 
which would lead to dynamical dark energy. 

If dynamical scalar fields are present, one naturally ex¬ 
pects them to couple to the rest of the model, unless a 
yet-unknown symmetry is postulated to suppress these 
couplings. In particular, a coupling to the electromag¬ 
netic sector will lead to spacetime variations of the fine- 
structure constant a —see 0 0 for recent reviews. There 
are some recent indications of such a variation 0, at the 
level of a few parts per million, which a dedicated Large 
Program at ESO’s Very Large Telescope (VLT) is aiming 
to test 0 0. 

Here, in the same spirit of EH , we discuss how as¬ 
trophysical and local tests of the stability of a can be 
used as additional tests of the underlying theories, in 
particular if one makes the ’minimal’ assumption that 
the same dynamical degree of freedom is responsible for 
the dark energy and the a variations. In this case ob¬ 
servational tests of the evolution of a directly constrain 
dark energy. The future impact of these methods as a 
dark energy probe has recently been studied in detail 
in E, [l3j].ffiarlier, more simplistic forecasts can also be 
found in [lj, Qr> .) Here we show how current data al¬ 
ready provides useful constraints. 

We start by reviewing the relation between a varying 
a and dynamical energy in the case (dubbed ’Class I 


* Carlos.Martins@astro.up.pt 
t Ana.Pinho@astro.up.pt 


models’ in 0) where both are due to the same dynamical 
degree of freedom. We will then focus on the case of a 
constant dark energy equation of state, wq, constraining 
this model with a combination of cosmological data and 
local and astrophysical tests of a. This choice is in the 
interest of simplicity, as it minimizes the number of free 
parameters and is frequently used as a fiducial model for 
forecasts; we leave a discussion of more general models 
for a subsequent, more detailed publication. 


II. VARYING a AND DARK ENERGY 

Dynamical scalar fields in an effective 4D field theory 
are naturally expected to couple to the rest of the the¬ 
ory, unless a (still unknown) symmetry suppresses this 
coupling E, El • We assume this to be the case for the 
dynamical degree of freedom responsible for the dark en¬ 
ergy. Specifically the coupling between the scalar field, 
</>, and the electromagnetic sector stems from a gauge 
kinetic function Bp((f>) 

C*f = -Ib f {<I>)F iiv F* v (1) 

which one can assume to be linear, 

B F {4>) = i - C«(0- M , ( 2 ) 

(with k 2 = 8irG): as pointed out in E the absence 
of such a term would require a <fi —> —(j) symmetry, but 
such a symmetry must be broken throughout most of the 
cosmological evolution. Local tests of the Equivalence 
Principle lead to the conservative constraint on the di¬ 
mensionless coupling parameter (see 0 for an overview) 

\(local\ < 10“ 3 , (3) 

while in El an independent few-percent constraint on 
this coupling was obtained using CMB and large-scale 
structure data in combination with direct measurements 
of the expansion rate of the universe. 
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With these assumptions one can explicitly relate the 
evolution of a to that of dark energy, as in M- The 
evolution of a can be written 

Aa^o-oo = _ (4) 

a ao 

and, since the evolution of the putative scalar field can 
be expressed in terms of the dark energy properties 
and w as 


1 + Wrj, 


30^ 


( 5 ) 


(with the prime denoting the derivative with respect to 
the logarithm of the scale factor), we finally obtain 

— (*)=C r Jmt{z){\+w+{z)) -^ 7 - (6) 

a Jo 1 ~r z 

The is assumes a canonical scalar field, but the argument 
can be repeated for phantom fields [18|, leading to 


— (Z) = -C r J3S1+(Z) ll+w+M-^ 
a Jo i- J- z 


( 7 ) 


the change of sign stems from the fact that one expects 
phantom filed to roll up the potential rather than down. 

In the present work we’ll focus on models with a con¬ 
stant equation of state wo, and will constrain them using 
the following datasets 


• Cosmological data: the Union2.1 Type la super¬ 
nova dataset fl9] and the compilation of Hubble pa¬ 
rameter measurements from Farooq & Ratra [201 ]. 


• Laboratory data: the atomic clock constraint on 
the current drift of a of Rosenband et al. 
which we can write as 


21J, 


= (-2.2 ±3.2) x 1(T 7 . (8) 

U o o. 


• Astrophysical data: we will use both spectroscopic 
measurements of a of Webb et al. [?] (a large 
dataset of 293 archival data measurements) and 
the smaller and more recent dataset of 11 dedi¬ 
cated measurements listed in Table [I] The latter 
include the early results of the UVES Large Pro¬ 
gram for Testing Fundamental Physics [$:, 9|, which 
is expected to be the one with a better control of 
possible systematics. 


III. OBSERVATIONAL CONSTRAINTS 

We now use the above datasets to constrain the dynam¬ 
ical dark energy model with a constant equation of state 
Wo ■ The behaviour of a will be determined by Eq. © for 
wq > —1 and Eq. 0 for w o < — 1. Our main interest is 


Object 

z 

A a/a 

Spectrograph 

Ref. 

3 sources 

1.08 

4.3 ±3.4 

HIRES 

[22] 

HS1549+1919 

1.14 

-7.5 ±5.5 

UVES/HIRES/HDS 

m 

HE0515—4414 

1.15 

-0.1 ±1.8 

UVES 

[23] 

HE0515—4414 

1.15 

0.5 ±2.4 

HARPS/UVES 

[241 

HS1549+1919 

1.34 

-0.7 ±6.6 

UVES/HIRES/HDS 

m 

HE0001—2340 

1.58 

-1.5 ±2.6 

UVES 

[251 

HE1104—1805A 

1.66 

-4.7 ±5.3 

HIRES 

[221 

HE2217—2818 

1.69 

1.3 ±2.6 

UVES 

m 

HS1946+7658 

1.74 

-7.9 ±6.2 

HIRES 

[22] 

HS1549+1919 

1.80 

-6.4 ± 7.2 

UVES/HIRES/HDS 

M 

Q1101—264 

1.84 

5.7 ±2.7 

UVES 

[23] 


TABLE I. Recent dedicated measurements of a. Listed are, 
respectively, the object along each line of sight, the redshift 
of the measurement, the measurement itself (in parts per mil¬ 
lion), the spectrograph, and the original reference. The first 
measurement is the weighted average from 8 absorbers in the 
redshift range 0.73 < z < 1.53 along the lines of sight of 
HE1104-1805A, HS1700+6416 and HS1946+7658, reported 
in [22] without the values for individual systems. The UVES, 
HARPS, HIRES and HDS spectrographs are respectively in 
the VLT, ESO 3.6m, Keck and Subaru telescopes. 


in obtaining constraints on the £—wq plane, and for this 
reason we will fix Hq = 70 km/s/Mpc and fl m o = 0.3 
(and assume a flat universe, so = 0.7). This choice 
of cosmological parameters is fully consistent with the su¬ 
pernova and Hubble parameter data we use. Moreover, 
we have explicitly verified that allowing Hq, fl m or the 
curvature parameter to vary (within observationally rea¬ 
sonable ranges) and marginalizing over these parameters 
does not significantly change our results: it is clear that 
the critical cosmological parameter here is wq itself. 

We therefore consider a 2D grid of ( and wo values, 
and use standard chi-square techniques to compare the 
models to the aforementioned datasets. Figure Q] shows 
the results of this comparison for the Webb et al. data 
(top panel), and for the Table Q] data (middle panel) in 
both cases, the constraints from the astrophysical data 
are shown by the thin red lines. The Webb data is not 
consistent with the null result [7jj, and we correspond¬ 
ingly find a one sigma preference for a non-zero coupling 
£ (with a negative sign for a canonical field, or a positive 
sign for a phantom held). However, the data is com¬ 
patible with the null result at two sigma. On the other 
hand, the Table [T] data is fully compatible with the null 
result. We note that in the former case the reduced chi- 
square of the best-fit model is XminWebb = 1-04, while 
in the latter case it is Xmin Table = 1-29; this may be an 
indication that some of the uncertainties in the Table |T| 
measurements are underestimated. 

For comparison we also show in the bottom panel of 
Fig. m in the same scale as before (and also in thin red 
lines), the local atomic clock constraint of Rosenband et 
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FIG. 1. One, two and three sigma constraints on the ( — wo 
plane from Webb et al. data (top panel), Table [I] data (mid¬ 
dle) and the atomic clock bound (bottom). In each panel the 
thin red lines correspond to the constraints from the astro- 
physical or clock data alone, the blue vertical ones correspond 
to the cosmological data (which constrain wo but are insensi¬ 
tive to () and the black thick lines correspond to the combined 
datasets. 



FIG. 2. One, two and three sigma constraints on the ( — wo 
plane from the full dataset considered in our analysis: Webb 
et al. data plus Table[I]data plus atomic clock bound plus cos¬ 
mological (Type la supernova and Hubble parameter) data. 
The reduced chi-square of the best fit is Xminjull = 0.97. 

al. [2ll | . For the models under consideration this trans¬ 
lates into 

1 Of / 

tt—= TCyS^^oll + w 0 | (9) 

£1Q OL v 

(with the — and + signs respectively corresponding to the 
canonical and phantom field cases), and it is clear from 
the plot that this is currently more constraining than the 
astrophysical measurements. 

The cosmological data we are considering is insensitive 
to (. (Strictly speaking, a varying a does affect the lu¬ 
minosity of Type la supernovas, but as recently shown in 
na for parts-per-million level a variations the effect is 
too small to have an impact on current datasets, and we 
therefore neglect it in the present analysis.) Naturally, 
the cosmological data does constrain wo, effectively pro¬ 
viding us with a prior on it. In Fig. [l] the cosmologi¬ 
cal data constraints are shown by the blue vertical lines, 
while the combined (cosmological plus astrophysical, or 
cosmological plus atomic clocks) constraints are shown 
by the thick black lines. Naturally, we can obtain tighter 
constraints by combining all datasets; this is straightfor¬ 
ward to do since the Webb et al. and Table [T] measure¬ 
ments are independent. The results of this analysis are 
shown in Fig. [2] 

Finally, in addition to constraints in the two- 
dimensional (-wo plane it is also interesting to obtain 
ID constraints on the coupling ( by marginalizing over 
the dark energy equation of state wq- The results of this 
analysis are shown in Fig. [3] We confirm that in the case 
of the Webb et al. dataset there is a one-sigma prefer¬ 
ence for a non-zero coupling, while in the other cases the 
null result provides the best fit. Significantly, the com¬ 
bination of all datasets allows us to obtain a non-trivial 
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FIG. 3. ID likelihood for £, marginalizing over wo . Plotted 
is the value of A\ 2 = y 2 ~ Xmim for cosmological + Webb 
data (blue dashed), cosmological + Table |T] data (blue dash- 
dotted), cosmological + atomic clock data (red dotted) and 
the combination of all datasets (black solid). 


constraint on At the two-sigma (95.4%) confidence 
level we find 

|C|<5xlO- 6 , (10) 

significantly improving upon previous constraints. As 
previously mentioned, the atomic clock measurement of 
Rosenband et al. currently provides tighter constraints 
than the astrophysical measurements. This new bound 
is the main result of our analysis. (Nevertheless, we also 
note that at three-sigma C, is unconstrained.) We can 
similarly obtain the ID likelihood for wq by marginaliz¬ 
ing over (. In this case we find at the three-sigma (99.7%) 
confidence level 

- 1.05 < w 0 < -0.94, (11) 

although this bound should be interpreted more cau¬ 
tiously given our assumptions on other cosmological pa¬ 
rameters. We leave a more systematic exploration of the 
relevant parameter space for a subsequent publication. 


IV. CONCLUSIONS AND OUTLOOK 

We have used a combination of astrophysical spec¬ 
troscopy and local laboratory tests of the stability of the 
fine-structure constant a, complemented by background 
cosmological datesets, to constrain the simplest class of 
dynamical dark energy models where the same degree of 
freedom is responsible for both the dark energy and a 
variation of a. In these models the redshift dependence 
of a depends both on a fundamental physics parame¬ 
ter (the dimensionless coupling £ of the scalar field to 
the electromagnetic sector) and background ’dark cos¬ 
mology’ parameters- for the simplest class of models we 
studied these are the dimensionless dark energy density 
fltf, and the dark energy equation of state wq. 

We obtained new, tighter constraints on the dimen¬ 
sionless coupling £ of the scalar field to the electromag¬ 
netic sector. We note that these constraints are cur¬ 
rently dominated by the atomic clock tests, which are 
only sensitive to the dark energy equation of state today. 
Thus a constant equation of state cosmological model is 
a reasonable assumption. Improvements in astrophysi¬ 
cal measurements will allow more generic models to be 
constrained. 

We have also pointed out how different currently avail¬ 
able astrophysical measurements of a (specifically the 
archival data of Webb et al. and the dedicated measure¬ 
ments of Tabled lead to somewhat different constraints). 
This highlights the importance of obtaining improved as¬ 
trophysical measurements of a (both in terms of statisti¬ 
cal uncertainty and in terms of control over possible sys- 
tematics), not only for their own sake but also because 
there can have a strong impact on dark energy studies. 
The next generation of high-resolution ultra-stable spec¬ 
trographs such as ESPRESSO and ELT-HIRES will be 
ideal for this task. A roadmap for these studies is out¬ 
lined in pi], and more detailed forecasts of the future 
impact of these measurements may be found in [13]. 
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